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DESCRIPTION 



[0001] The object of this invention is a process for reconstruction of a tomographic image, 
of the analytical type, and in which improved modelling of the movement of the object is 
resorted to so as to reduce the artefacts of the image. 

10 [0002] The reconstruction of images in tomography implies the utilisation of radiation 
passing through the object. In numerous cases, the radiation is irradiation radiation partially 
attenuated by the points of the object it passes through; it also happens that the radiation 
lines are artificial and correspond simply to collimation lines of the detectors, which register 
an emission of particles in the object. The first of these processes are radiography processes, 

15 and the others are emissive processes, reserved more for living beings and in which the 
radiation is produced by an emissive body previously absorbed. Although the two families 
of processes are completely different, the majority of the processes for reconstruction of 
images apply to both, and this is also the case with the invention. 

[0003] Divergent radiation is often preferred for more easily enclosing the object studied. 

20 In the case of irradiation radiation attenuated by the object, this is achieved by having a 

source emissive at a punctiform focal point near the object and a network of detectors on the 
opposite side of the object, all of which are collimated towards the source; in the case of 
radiation emitted by the object, a similar network of detectors is collimated to a un focal 
point which corresponds geometrically to the preceding punctiform source, but to any 

25 material object. In all cases, the focal point and the network of detectors are moved about the 
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object by taking successive views thereof, and for each of the views, the detectors measure 
totals of the attenuation or emission property along collimation lines, known as projections 
of the image of the object. When a sufficient number of views has been taken, there is a 
large number of projections criss-crossing through the object. The mathematical problem 
5 known as system inversion provides the attenuation or emission property which serves to 
form the image at each of the points of the object from sums of properties on the projections. 
It can be represented by a linear system of equations where the values known at the outset 
are the measurements taken by the detectors at different locations and the unknown ones are 
the values of the property at different points of the object. Certain methods of reconstruction 

10 known as algebraic effectively resort to inversion of this system; yet there are other methods 
known as analytical, where the value of the property at each point is calculated directly from 
a mathematical combination of the projections. The present invention is one such. 
[0004] An important domain of tomography is the study of living beings for medicine. As 
several successive views must be taken, the object in question is capable of moving between 

15 them according to a general law of evolution. The different views thus represent different 
states of the object, and artefacts of reconstruction can appear only on the reconstructed 
image. 

[0005] Processes for avoiding this permanent problem consist of taking all the views at 
once with as many sources and networks of detectors, or taking views of the object only in 
20 states identical to the latter, which is possible if its movement is periodical as are certain 

physiological activities such as heartbeat or breathing; the first of these processes is however 
expensive and the second is delicate to implement conveniently. 
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[0006] This is also why numerical models of the movement of the object reflecting the 
general law of evolution of the object were incorporated in the process of reconstruction. 
The model generally goes back to a collection of shift fields of the points of the object at 
respective times where the views are taken. This can be established by other measures or by 
5 certain hypotheses. Several genres of models have already been put forward, but they are 
generally insufficient. US patent A 5 287 276, applied to the study of the thoracic cage, 
considered translations and dilations of its content, but this model is insufficient for other 
movements, such as that of the heart, which comprises torsion. In addition, the sets of 
projections utilised to apply the inversion formula are not very well known, which can allow 
10 new artefacts. 

[0007] The process of the invention offers the advantage of using a more elaborate but also 
very simple model of the movement of the object, comprising new classes of movements, 
and returning to a new and perfectly exact inversion formula. 

[0008] In addition, it is easy to obtain the sets of projections necessary and only necessary 

15 for reconstruction of the points of the image. 

[0009] Generalisations to more complex situations are also proposed. 
[0010] In its general form the invention concerns a process for reconstruction of a 
tomographic image of an especially mobile and deformable object, the image being a set of 
values of a property taken by points of the object, comprising the use of: divergent radiation 

20 from a focal point and passing through the object, the focal point being mobile about the 
object; an analytical model of mobility and deformation of the object defined for each 
position of the focal point; and an analytical calculation process for obtaining said values 
from totals of the values of the property along projection lines leading to the focal point and 
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passing respectively by the points; characterised in that the model is improved, and is a 
variable combination being acquired, this combination comprising translations, rotations and 
homotheties of the object from an origin, and in that the process of analytical calculation 
comprises the following stages: 
5 - weighting of the measurements, this weighting being dependent on the analytical model 
of mobility and deformation of the object; 
- derivation of the measurements weighted following the trajectory of the focal point 
considering a direction adapted to the model, this direction being kept constant, and 
obtaining modified measurements; 
10 - retroprojection of the modified measurements. 

[0011] Prior to weighting, filtering of the measurements acquired by a Hilbert filter is 
often applied, since it is usual in other processes. The weighting and derivation then apply to 
the filtered measurements. 

[0012] The following figures are introduced to explain the invention: 
15 Figure 1 illustrates taking the measurements of the object in question, 

Figure 2 illustrates an equivalent diagram, the object being maintained in a state of 
reference, 

Figures 3a, 3b and 3c illustrate a complex case of deformation of the object, 
Figure 4 is similar to Figure 2 in the complex case, 
20 Figure 5 illustrates taking measurements with conical radiation, and 

Figures 6, 7 and 8 are organigrams of three embodiments of the invention. 

[0013] Reference is now made to the figures. The first of them illustrates the classic 
situation in tomography of a source S of radiation and a tuned mobile network of detectors D 
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on circular trajectories concentric to opposite positions around an object E to be studied. 
Radii R join the source S to the respective detectors of the network D. Under consideration 
here is a problem of geometric plan, which corresponds to the so-called fan-shaped parallel 
measurement conditions, where study of the object E is made by superposed slices s. The 
5 processes implying conical radiation are also common and will be examined hereinbelow. 
The invention was designed to divergent radiation in general, which include the latter also. 
[0014] The trajectory T of the source S and opening of the radiation bundle are generally 
selected such that the latter includes the whole section of the object. 
[0015] We should consider a point P of the object E. A sole radius originating from the 

10 source S passes through it when a view is taken, and other radii R2, R3, etc. likewise pass 
through it for other positions S2, S3, etc. from the source S when other views are taken. 
These beams passing through the point P are utilised to determine the image of the point P in 
the inversion process. It is useful that the radii being considered are oriented in directions as 
varied as possible, so as to describe all the straight directions possible for each point of the 

15 object E. 

[0016] It should now be considered that the object E is moved and deformed. The point P 
moves inside the trajectory T and the radii R2, R3, etc. to be employed pass through 
respective positions P2, P3, etc. distinct from the point origin P. The problem of 
reconstruction to be resolved is identical to the artificial configuration of Figure 2, or the 
20 point P was taken as immobile reference and where everything happens as if the source S 
followed a trajectory T having an irregular form and to the numerical definition noted by 

OS =T A (5(A)) (and 5(A) in the real configuration of Figure 1). It is obvious that this 
artificial diagram of the problem aids in its comprehension. According to the invention, the 
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model of displacement and deformation selected for the object E is improved, composed for 
example from a variable combination over the course of acquisition of translations, rotations 
and homotheties, according to the formula (1) 



*o= r A(*) = 4i*+ 5 A = 



a u (A)a n (X) 



x+ 



b 2 (X) 



(i) 



a 2X (X)a 22 (X) _ 

5 where x 0 is the vectorial position of the point P relative to a reference such as the point O at 

the reference instant selected to carry out the calculations and reconstruction of the image, x 
is the position of the 

point 2 at another instant and especially a view-taking instant, and the coefficients a and b 
dependent on the time. X is a general parameter of the process, which varies between a 
10 minimal value X min and maximal value X max; the time is fixed by a monotone function de 
X, (optionally by X itself). The reference instant where the image is reconstructed 
corresponds to X = 0. The attenuation or emission property to be calculated is noted f, and 

for the point P fixed by the vector x it is noted fX{x) at any X instant and f 0 (x 0 ) at the 

reference instant, where x is noted as x 0 . 
1 5 [0017] The value of projection of the property measured over the entire radius R going 
from the source to a detector of the network D via the object E is given by the formula (2), 

where a is a unitary vector, and this projection is noted as a function of this unitary vector 
a and of the parameter X: 

20 gx&,S)=l R dtfA*W + t*) ( 2 ) 

where t is a parameter expressing the progression of the projection on the radius 

and taking the values included in the set of positive real numbers R. 



0 
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[0018] It is usual to apply a Hilbert filter to the projections in analytical methods, which is 
done here also, and the filtered projections are noted gHX according to the formula (3) 



5 1 

[0019] Where n is any vector, S 1 is the sphere unit in tridimensional, Iih is indicated by the 
5 formula (4) 



[0020] To create the formulas the specialist could use either a classic operation of discrete 
convolation which will introduce a discrete version of the Hilbert filter, or a multiplication 
operation in the domain of Fourier, which will introduce an apodised version of the Fourier 
10 transform of the Hilbert filter. 

[0021] The inventors have also established the following equation (5) 



projection straight line at the instant X passing through the mobile point corresponding to the 

15 point P, n is still any vector and or P 0 h designates parallel projections ¥ 0 h ( n , s), obtained 
by rearranging the projections effectively measured on the object in the reference state and 
filtered by the Hilbert filter according to the following formula (6) 




h H (s)=- \daisign{G)e lixas . (4) 



-OO 



|det^| 

where the matrix-vector product Aln is the direction orthogonal to the 



Po H {n,s)= [ds'P 0 (n,s')h H (s-s'), 



s being a parameter giving the distance at the origin. 



20 The inversion formula carried out is thus the following formula (7) 
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where « * designates a unitary vector orthogonal to a radius linking the point P 
in question to the source S in the illustration of Figure 2, and where g F x is given in the 
following formula (8) 



[0022] The specialist could use a method of difference finished on two or three points to 
create this derivation step. 

[0023] This derivation is performed according to the parameter X of the trajectory of the 
focal point (source S). It is specifically adapted to the improved deformation and mobility 

10 model as it applies to a direction A T p n by considering that the direction n is kept constant. 

So it is not the direction orthogonal to an acquired radius, but the direction orthogonal to the 
equivalent radius of the artificial geometry, which is to be kept constant. 

[0024] A T p n is the direction orthogonal to the projection straight at the instant X passing 

through the mobile point corresponding to point P. 
15 [0025] This inversion formula takes deformations of the object E into consideration and 
comprises, relative to other formulas, set in ordinary cases, a weighting as a function of the 
deformation of the object (by the determinant of AX) and of the position of the trajectory (by 

the standard between x 0 and T/l(a(/l)).These barely constrained conditions are explained by 
an improved space of transformation, the projection straight lines of Figure 1 remain straight 
20 lines in the artificial geometry of Figure 2, such that the numerical problem can be resolved 
analytically. 
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In formula (7) 



||* 0 -r x (5(A))|| 



the limits of the integral A r (x 0 ) designate a minimal set of positions X of the source S, such 
that on the object E in the reference state the directions of the straight lines linking x 0 to 
5 r>i(tf (>!)) cover the entire interval of a semi-turn of trajectory without redundancy. If care is 




(7) could be classically made discrete by the specialist by a formula of trapezes. 

[0026] It is shown that the improved deformations do not sufficiently model all the objects 

which have to be studied in practice. Figures 3a, 3b, 3c, and 4 are offered as a revision of the 
10 reasoning of figures 1 and 2. The object E evolves in complex fashion between the positions 

indicated successively in states XI, X2 and X3 or the position of the source S was also 

indicated along with the radius Rl, R2 or R3 leading to point P. 

[0027] In returning the object E to a reference state, Figure 4 shows that the radii 

designated R'l , R'2 and R f 3 leading to point P and corresponding to radii Rl , R2 and R3 
15 (passing through the same points of the object E) are no longer rectilinear. The formula (7) 



is thus no longer directly acceptable. However, it is approached by approximations, 
20 supposing improved deformation of the object E particular to each point P and valid around 
it. The deformation matrix corresponding to A\, of the formula (1) 




taken to thus limit the integral, the reconstruction formula is perfect. The integral of formula 
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a n (X)a X2 (X) 
a 2 ,(X)a 22 (A)^ 



x+ 



b 2 {X) 



(1) 



is thus calculated by the formula (9) 



^(r-'(*o))= 



9 (r^.r-foi-^tf.r^fo)) 



cbt. 



cbt, 



5 where T xl (A,x) and T x2 (A 9 x) are the compounds according to the main directions xi and x 2 
of . B\ does not need to be calculated. The estimated shifts of all the points of the 

object E are thus the coefficients T(X 9 x) which can concretely be introduced to shift cards 
which are read and executed at the moment of inversion. 



10 [0028] At each point P, it will still be necessary to determine the set of projections 

A r (jc 0 ) which will be used for inversion, according to the abovementioned principle that an 
angular interval of a semi-turn must be covered by these projections. As the radii of the 
artificial geometry R f l, R'2 and R'3, etc. are no longer rectilinear, the directions of their 
tangent to the point P of intersection are considered. 

15 [0029] The inversion formula can thus be written according to the formula (10) 



where gpx is expressed by the formula (11) 
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[0030] In addition, n\T' ] [xo U lis the direction orthogonal to the straight line acquired at 



5 [0031] Finally, L(x 0 -T a (a(A))is the distance from point P to the image of the source on 
the reference object, along the virtual radius R\ 

[0032] To the detriment of quality, but in the interests of greater simplicity, the term 
L(x 0 -T A {a(X)) can be replaced by \\x 0 -T x (a(X))l 



10 [0033] The improved model can be defined according to the real deformation of the object 
by an approximation according to an approximation criterion such as the criterion of lesser 
squares, the criterion of minimisation of the L 1 and L 2 standard, optionally complete by 
regularising on the gradient or the Laplacian. 

[0034] To date interest has been shown in reconstructions of the object E by slices, in fan- 
15 shaped parallel conditions. The preceding process can be extended with acquisitions by 
conical radiation under the conventional conditions of Figure 5, or the object E is 
represented in its three-dimensional entirety, also as for the network of detectors D, which is 
composed of a series of layers of detectors, each similar to that of the preceding figures; a 
radius R can be expressed by three parameters, namely the position X of the source S, an 
20 angle y which makes the radius R in the plane of the trajectory T relative to the central axis 
X of the bundle, and a cote q marking the layer of detectors in which the radius terminates. 



the instant X passing through 



T l (x 0 ) ; /?(«(r l (jf 0 ),/l))is a factor associated with 



deformation of the lengths on the radii, according to the formula (12) 
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The angle 7 replaces the direction vector a previously used in the formulas out of 



convenience. 



The following formula (13) 



L(x 0 -rM*)J 



;8fx 



5 gives the inversion which is undertaken in these particular geometric conditions and by 
making a Feldkamp approximation commonly used for treating the conical projections. The 
measurements are in this case multiplied by cos A, which is classic weighting compensating 
the oblique character of the acquired radius. This angle A is given in Figure 5 and illustrates 

the angle of the radius R with the plane of the trajectory T. In addition, l[x 0 -Tx{a{jijjj is 




10 the distance from x to the image of the source S on the object E in the reference state along 



direction orthogonal to the straight line acquired at the instant X passing through Y l [x0) on 
this plane. Finally, it is apparent that q[X 9 T~ l (x 0 )) is the side on the network of detectors D 
of the straight line, acquired at the instant X, passing through the mobile point corresponding 
15 to;c0. 

[0035] A particular process, first described in French patent 01 0791 8 filed under number 
can still be applied here by way of some transformations. This process applies to rapid 
acquisitions made around the object by making several turns and consists of a reconstruction 
of the animated object by evaluating its contents by blocks of projections taken on a fraction 
20 of a turn only of the source. The inversions of projections made on the blocks produce sub- 
images which are incorrect since they comprise only part of the measurements, but are 
obtained very rapidly and with the same rapidity create the law of displacement or 




the virtual radius R', in projection on the plan of the trajectory T, and n\T l (x6),A/j is the 
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deformation of the object by comparing homologous sub-images, taken for the same 
positions of spaced sources of a complete turn or a semi-turn; the sub-images are 
reconstructed at a reference instant for each of the groups of blocks and finally combined to 
give the complete image of the object. Here, the limits of the blocks are given in virtual 
geometry. The process is not otherwise modified. 

[0036] Finally, it is possible to improve the process to study an object subjected to 
periodical phenomena. Only those which are taken at phases of the phenomenon similar to 
those of the reference instant are selected as study blocks. 

[0037] Reference can be made to the organigrams of Figures 6, 7 and 8 to understand the 
invention. They detail three embodiments of the invention, or reconstruction is done 
respectively with compensation of movement; compensation of movement and time 
compensation; and compensation of movement, time compensation and consideration of 
periodicity of the evolution of the object. What has been previously described is applied in 
these processes. The voxels are evidently the points of the image of the object considered in 
the reconstructed image. 
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